Relation between Zitterbewegung and the charge conductivity, Berry curvature and 

the Chern number of multi band systems 
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We show that the charge conductivity for impurity free multi band electronic systems can be 
expressed in terms of the diagonal and non-diagonal elements of the Zitterbewegung amplitudes 
while the Berry curvature and the Chern number is related only to the diagonal elements. Thus, the 
phenomenon of the Zitterbewegung can no longer be viewed just as an interesting consequence of 
quantum physics but it has also an experimental relevance. Moreover, through several examples we 
demonstrate how efficient our approach is in the analytical calculation of the charge conductivity. 
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Introduction. — The Zitterbewegung predicted origi- 
nally by Schrodinger for Dirac electron is a 'trembling' 
or in other words a rapid oscillatory motion of the center 
of the free wave packet for relativistic electron Most 
recently, Schliemann et al. [2\ predicted the Zitterbewe- 
gung in spintronic systems where the experimental ob- 
servation of the effect is more realistic due to the much 
smaller frequency of the oscillatory motion. Thus, the 
Zitterbewegung can, in principle, be observed not only 
in the relativistic regime but for spintronic systems [2J 
and graphene [3| as well. This seminal paper has initi- 
ated many other works with an aim to demonstrate the 
appearance of the Zitterbewegung not only for Dirac elec- 
trons but for quasi-particles in condensed matter physics 
(see references in our recent work on the general theory 
of Zitterbewegung where alternative expressions for 
the Zitterbewegung amplitudes and explicit forms of the 
position operator for several systems are presented). 

In connection with graphene Katsnelson has already 
pointed out that the Zitterbewegung resulting in an os- 
cillating term in the current operator is responsible for 
the nontrivial behavior of the conductivity at zero tem- 
perature and zero chemical potential [5] . In this paper we 
show generally that the charge conductivity for a impu- 
rity free multi band system is related to Zitterbewegung. 
Such a relation is expected since in the Kubo formula 
for the charge conductivity one needs to calculate the ve- 
locity operator in Heisenberg picture which includes an 
oscillatory term due to the Zitterbewegung in the case of 
multi band systems. Moreover, we clarify the relation- 
ship between Zitterbewegung and the Hall conductivity 
of insulators with non- interacting Bloch electrons [H, @] ■ 

To show the subtle relation between the Zitterbewe- 
gung and the charge conductivity we start with a multi 
band system described by the most general matrix Hamil- 
tonian in a Bloch wavefunction basis: H a b(k), where 
a,b = 1,2, ■■■N are the band indices (here N is the 
number of bands of the system). Here each matrix ele- 
ment H a b(k.) is a diffcrentiable function of the wave num- 
ber k corresponding to the Bloch states. In Ref. [1] we 



calculated the time dependence of the position operator 



x(t) 
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of the quasi-particle by decom- 



posing the Hamiltonian into a sum of projection opera- 
tors: H = J2 a £-a(k)Q Q (k), where E a (k) are the distinct 
eigenvalues of the Hamiltonian at a given wave number 
k, and Q a (k) are projection operators (N x N matrix 
operator) satisfying the usual relations: Q a Qb = 5 a b Qa 
and ^2 a Q a = In , where In is the N x N unit matrix. 
Note that a = 1, 2, • • • , s, where s < N (for degenerate 
case s < N) . We found that the time dependence of the 
position operator x(t) in Heisenberg picture becomes 

x(i) = x(0)+Wt + J2 zab (e lul « bt - 1) , (la) 

a, 6 
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W = ft M- Qa > Z =t Qa — , (lb) 



9k 



dk 



and uj„ 



E \ Eh . Here W is the drift velocity and Z ab 
are the Zitterbewegung amplitudes. This is a general re- 
sult for the phenomenon of the Zitterbewegung (a simple 
derivation of the above result and explicit examples are 
presented in our recent work Q). In what follows, these 
Zitterbewegung amplitudes Z ab will play a crucial role in 
the charge conductivity and the Chern number. 

Using the Kubo formula we show that the frequency 
dependent charge conductivity (often called optical con- 
ductivity) of a multi band system can be expressed in 
terms of the Zitterbewegung amplitudes Z 
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n F (E a - fj) - n F (E b - fj) 
iv m + E a - E b 



(2a) 
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(2c) 



Here the Zitterbewegung amplitudes Zf' (i denotes the 
components x,y,z) are given by Eq. (jTbJ) , v m — 2mrn/f3 
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(m is an integer, (3 = l/(k-BT)) are the bosonic Mat- 
subara's frequencies, n F (E) = l/(e^ E + 1) is the usual 
Fermi distribution, /i is the Fermi energy, Im{-} is the 
imaginary part of the argument, f stands for the conju- 
gate transpose and finally V is the volume of the sam- 
ple. The trace is taken over the band indices. Note that 
K a b{z) — Kba(-z) and K aa (z) = for z ^ 0. Equation 
(|5J) is one of our central results in this paper. Here we 
omit the calculation of the Drude peak, and focus on the 
inter-band contribution. One can see clearly from Eq. ([2]) 
that the Zitterbewegung amplitudes describing the inter- 
band interference manifest in the inter-band contribution 
of the charge conductivity. The Zitterbewegung ampli- 
tudes are indeed probed every time when one measures 
the conductivity of a scattering free multi band system. 

Derivation of Eq. (0). — First consider the response 
function for operators A and B given by 



V Jo 



(f A{t) B(Q)) e w ™ T air 



where G(k, z) = [z + fx — is the one-particle 

Green's function, T is the time ordering operator and 
LJ n = (2n + l)7r//3 are the fermionic Matsubara's frequen- 
cies (n is an integer). Using the projector decomposition 
of the Hamiltonian H = EaQa the Green's function 
takes the following form: G(k, z) = z+^-e (k) • Then 
substituting this Green's function into Eq. (|3]) one can 
find 

6abK) - ~ £ Kbaiivm) Tr( A Q a BQ b ), (4a) 

k 

K ah {iv m ) =pJ2j 
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iui n + iv m + fi- Eb iu„ + H - E a 
n F (E a - (j,) - n F (E b - fx) 



E a — Eb 



(4b) 



where in the last step of the calculation of the function 
Kab{iv m ) we have used the usual summation technics 
over the Matsubara's frequencies Q. 

Now using Eq. Q the current-current correlation func- 
tion with current operator J = ^ (in units of e/h which 
is taken into account in the expression of the conductiv- 
ity) reads 



Uij(iu m ) 



7 < 5 » 



k a,b 



Now using the well-known relations Qa^Qb = 
Sab ^Qa + (E b - E a ) Q a ^ and Q a Q b = S ab Q a , and 
the fact that K aa (z) — for z =^ 0, it is easy to obtain the 



current-current response function Ily given by Eq. (|2bl) . 
Finally, the conductivity ((2a)) can be calculated by ana- 
lytic continuation iv m — > fvjj + iS, where 6 is a positive 
infinitesimal. 

The Berry curvature and the first Chern number. — 
Now, we show that the first Chern number characterizing 
the Hall conductivity of insulators with non-interacting 
Bloch electrons [f| 0] can also be expressed in terms of 



the Zitterbewegung amplitudes (lb). 

From Eqs. <j5j) and ([2a)) we find that the intrinsic Hall 
conductivity (here we focus on two dimensional systems) 
for uj — > (dc conductivity, for j ^ I) is 



<Tjl(L> = 0) 



nS»>(k) 



h A 



£J>n<7>(k) 



k a 



2 ^E 

b^a 



b - c.c. 





{E a - 


-E b f 




\dQ a 






\Qa 


_dk J 


dh 


) 



(6a) 

(6b) 
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where n a = n F (K(k) - Li), ^f(k) is the Berry curva- 
ture (Eq. (|6b|) can be casted to that derived by Thou- 
less et at. in Ref. [Uj]), A is the area of the sample 
and c.c. and [•, •] stand for the complex conjugation 
and the commutator, respectively. To get (|6c]l we used 
again Q a ^Q b = S ab ^Q a + (E b - E a )Q a ^ and 
Qa + J2b^a Qb = In- Now replacing Q a by Q a = Q 2 a in 
c]) we obtain a very simple expression for the Berry cur- 



vature in terms of the Zitterbewegung amplitudes: 



n^(k) = 27riTr 



(7) 



It is interesting to note that in contrast to the Zitterbewe- 
gung, where only the non-diagonal elements Z ah (a ^ b) 
appear in (|lb|) . in the Hall effect only the diagonal ones 
Z aa play the role. However, in the charge conductivity 
for finite frequencies given by Eq. the non-diagonal 
elements Z ab are present. 

The Hall conductivity for band insulator in which the 
Fermi energy \i is located inside the energy gap be- 
tween conduction and valence sub-bands and at zero 
temperature reads as * B = ^ where C x = 

\/(2A) EkEE^^^'^jfC 4 ) i s the nrst Chern number 
(here Sji is the fully antisymmetric tensor and the sum- 
mation is assumed on indices j and I) and it can also be 
expressed with the Zitterbewegung amplitudes as 



Cl = -^-/^£^Tr( 

« a 



Zf a Z l 



aaf 



(8) 



Note that for the first Chern number one can ob- 
tained the same result starting from Eq. (3) given by 
Avron et at. in Ref. 12 1. 
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An important mathematical theorem. — Now, we recall 
a less known mathematical theorem which enables us to 
calculate the projector operators Q a = \a)(a\ without 
calculating the eigenvectors \a) of the Hamiltonian H. 
Let H be an N x N hermitian matrix with s < N dis- 
tinct eigenvalues, E a , . . . , E s . Then the matrix H can 
be decomposed in terms of projector matrices as H = 
J2 a E a Qa, where the projector matrix Q a for a — 1, . . . , s 
(in the mathematical literature called Frobenius covari- 
ant [3]) is given by 



where k = (k x ,k y ,k z ) is the wave number and S = 
(S x , S y , S z ) represents the spin operator with spin 3/2, 
while 7 12 and m are parameters of the model. The 
Hamiltonian can be expressed in terms of the projection 
operators Q+ and Q- as Q 

H = £_(k)Q_(k)+£: + (k)Q + (k), (13a) 

Q + (k) = \h- -^(kS) 2 , (13b) 
o 2 k 

Q_(k) - / 4 -Q+(k), (13c) 



Qa 



n 

6=1 



1 



Ea-E b 



(H-E b I N ). 



(9) 



The proof of this theorem is based on the Cayley- 
Hamilton theorem 



13|,ll4|. 



Applications. — In the following we show a few exam- 
ples how the phenomenon of the Zitterbewegung is re- 
lated to the charge conductivity or the Chern number 
for specific multi band systems. 

i) Consider the most general non-interacting two-band 
model 0, [i| with Hamiltonian: 



H = e(k) h + h(k) er, 



(10) 



where I2 is the two by two unit matrix in spin or pseudo 
spin space, and the system is characterized by the one- 
particle energy dispersion e(k) and an effective mag- 
netic field h depending on the wave number k, while 
c = (oxj °~y,o~z) is a vector formed from the Pauli matri- 
ces. The two eigenvalues are E±(k) = e(k) ± h, where 
h = V 7 !?, and the two corresponding projectors obtained 
from Eq. (|9]) arc Q± = | ( I2 ± h tr\ where h = h/h is a 
unit vector. The conductivity for impurity free samples 
can be obtained from Eq. ^ and after a little algebra 
the response function becomes 



Hij(iv m )-- 



£7- 



I dh a dhp ~ 



where n± — np(E± — fj,) and we sum on any repeated 
index. For two dimensional samples this result agrees 
with that obtained by Bernevig [9(. 

ii) For a two-band model with Hamiltonian (|10p us- 
ing Eq. ([5]) we easily find the first Chern number: C\ = 

~ Tk J d 2 k h • (^-§jr- x -jj^-j, which is a well-known re- 
sult 0,0. 

iii) We now consider the Luttinger-type systems [(| [l5[ 
for which the Hamiltonian is given by 



2m 



7i + 2 



k 2 - 2 72 (kS) z 



(12) 



where 7 4 is the 4x4 unit matrix, and the double degener- 
ate eigenvalues are E±(k) = 71 ^ 2-72 (Kk) 2 corresponding 
to the light-hole (+) and the heavy-hole (— ) bands. The 
projection operators Q± can be obtained from Eq. ([9]). 
In earlier calculations of the response function the 
SO(5) Clifford algebra has been invoked [1, In our 
approach the response function can be obtained without 
using the Clifford algebra. Indeed, it is easy to calcu- 
late the current-current response function using Eqs. (jlb[) 
and (I2bp or Eq. ([5]), and the commutation relations 
[Sj, Sf.] — iejkihSi for the spin operator S. The Zit- 
terbewegung amplitudes have already been calculated in 
Ref . Q . After some algebra we have 



TLij(ivm)= 



24h 6 ~f. 



, (n + — Ti-) (k 2 Sa — kikA 
d 3 k— ' y n 3 ^-,(14) 



d J k 
(2^F 



The 



where i,j — x,y,z and we used ^ k — > V J 
integration over the polar angles of k can be done analyt- 
ically and at zero temperature we have the same result 
as that, e.g., in Ref. [9j. 

iv) Consider the spin-orbit interaction in two- 
dimensional electron gas in a fully symmetric quantum 
well investigated recently by Bernardes et al. [16j. The 
Hamiltonian of this system is given by 



1 2m 



H 



V 



+ e e 
ii]k + 





-irjk- 
2 

- + 








irjk + 






-ink. 



\ 



+ £e / 



(15) 

where k± = k x ± ik y , while e ej0 , 77 and m are parameters 
of the model. The two double degenerate eigenvalues 
ar e E± = e k ± f k , where e k = h 2 k 2 /(2m) + e+, f k = 

\J £ 2 - + k 2 ?7 2 , e± = (e e ±e )/2. Again one can show that 
H = E + Q + + where the two projection operators 

obtained from Eq. © arc 



( ^ink- 
±ink + h k 

\ ±ink+ 











^fink- 



\ 



(16) 
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and h k = ±e 



+ fk- To calculate the current-current 
response function we again use Eqs. (jlbl) and ([2]). After 
some algebra we have 



2 f d i 

n«(«^)=2,y — 



d 2 k (n+ - n-) {fl5 lJ - kikjrj 2 ) 



fk 



ft 



,(17) 



where i,j = x,y. Performing the integration over the 
polar angle of k one can see that IIy(ii/ m ) is a diagonal 
matrix. To our best knowledge this result is new in the 
literature. The charge conductivity can be obtained from 
Eq. (j2"aj) and it will be published elsewhere. 

v) For single layer graphene the most general Haniil- 
tonian in tight-binding approximation can be given by 
a 2 by 2 matrix in which the diagonal elements Haa 
and Hbb include second, forth, etc. neighbor hopping 
terms, while the off-diagonal elements Hab and Hba 
contain the first, third, etc. hopping terms [13] • Thus, 
the Hamiltonian can be mapped to Eq. (fTU|) . Note that 
the same is true even for a strained graphene. The re- 
sponse function can be obtained from Eq. (JTTJ). In par- 
ticular, taking into account only first-nearest neighbors, 
we have H A a = H B b = £o and H A b = H* BA = /(k), 



where /(k) = 70 (1 + e~ lkai + e~ lka2 j, and 70 and e are 
parameters of the model, and ai and a2 are the unit vec- 
tors of the unit cell in the honeycomb lattice [TJj]. The two 
eigenvalues of the Hamiltonian are E± = Eq ± |/(k)| and 
h = (Re{/(k)},-Im{/(k)},0)/|/(k)|. Using Eq. (m 
the current-current response function is 



A V l/(k)| 2 



Ftf(k),(18a) 



Mf|}Mf|} 

l/(k)| 



(18b) 



where the summation in k is over the entire Brillouin 
zone of the honeycomb lattice, A is the area of the sam- 
ple, Re{-} is the real part of the argument and * denotes 
the complex conjugation. Note that Eq. (fT5|) is valid 
not only in the usual Dirac cone approximation. Our re- 
sult agrees with that obtained by Zhang et al. [l9[ and 
by Yuan et al. [20]. However, according to our numeri- 
cal calculations it slightly differs from that obtained by 
Stauber et al. [2l| for the high frequency region. 

vi) Finally, we give at least one example in which the 
projector decomposition of the Hamiltonian involves not 
only two projector operators (as in the previous cases) 
but four projector operators. Such a system is, e.g., 
the bilayer graphene 22, [23J. Our general framework 



presented in this work for calculating the charge con- 
ductivity can also be applied to bilayer graphene. The 
Zitterbewegung amplitudes have already been given in 
Ref. 0] for bilayer graphene excluding trigonal warping. 
Using this result we obtained the same analytical expres- 
sion for the frequency dependent optical conductivity as 



that by Nicol and Carbotte in Ref. 24] using the spec- 
tral function representation of the Green's function. We 
would like to stress that our approach for calculating the 
optical conductivity is a convenient and very efficient 
method even for more complex systems. The study of 
strained bilayer graphene with/ without trigonal warping 
is in progress. 

Conclusions. — In this paper we derived an explicit ex- 
pression for the charge conductivity, the Berry curvature 
and the Chern number in terms of the Zitterbewegung 
amplitudes. Our results show that the Zitterbewegung is 
not just an interesting phenomenon in quantum physics 
but it is closely related to measurable effect. 
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